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We show that a composite-field (diatom) Goldstone state is expected in a dilute Bose gas for 
temperatures between the Bose gas critical temperature where the atom Bose-Einstein condensate 
appears and the temperature where superfluidity sets in. The presence of superfluidity is tied to 
the existence of a U(l) charge-two diatom condensate in the system. By promoting the global U(l) 
symmetry of the theory to a gauge symmetry, we find that the mass of the gauge particle generated 
through the Anderson-Ifiggs mechanism is related to the superfluid density via the Meissner effect 
and the superfluid density is related to the square of the anomalous density in the Bose system. 
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Ultracold atomic gases provide a fertile testbed for 
many-body theory. Tabletop experiments can be de- 
signed to take advantage of the ability of tuning the 
strength of the inter-particle interactions and provide 
new insights in the role played by correlations in diverse 
physics problems such as high-Tf, superconductors and 
the equation of state for neutron matter. To this growing 
list of applications, we would like to add the Higgs mech- 
anism, which is of great interest because of the ongoing 
experiments at the Large Hadron Collider and their po- 
tential implications to the physics of the Standard Model 
of elementary particles. In addition, the Higgs mecha- 
nism in ultracold atomic gases is particularly intriguing 
because recent advent of synthetic magnetic fields allow 
for the study of electromagnetic-like effects in systems of 
neutral atoms [1]. 

Recently we reformulated the theory of dilute Bose 
gases in terms of normal and anomalous auxiliary-field 
(AF) densities [2] . In the leading order of the AF loop 
expansion (LOAF) we find two critical temperatures in 
the phase diagram: the critical temperature Tc, where 
the atom Bose-Einstein condensate (BEC) appears first, 
and a temperature T* > Tc, which indicates the onset of 
superfluidity and the appearance of a diatom condensate 
in the system. In this paper we will show that between 
Tc and T* the system supports zero-energy and zero- 
momentum excitations, which correspond to composite- 
field (diatom) Goldstone states. Introducing a formal 
U(l) gauge vector meson into the system, we find that 
the mass of the gauge particle generated through the 
Anderson-Higgs mechanism [3] can be related to the su- 
perfluid density via the Meissner effect. 

In the LOAF theory of dilute Bose gases [2, 4] the 
classical action is given by = J d^a;£[$], where 
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^[^]-2a[ 



x\x)-\A{x)\']-V2x{x)\Hx)\' (1) 
[A*ixmx)r+A{xWix)r] 

■ [ (j)* (x) h (t)ix) + (j){x) h* (t>* {x) ] , 



with h = ihdt + 7V^ + fi and 7 — h"^ /{2m). The cou- 
pling constant. A, is related to ao, the s-wave scatter- 
ing wave [4]. Here, x{x) and A{x) are real and com- 
plex fields related to |^(a:)p and 4>^{x), respectively, and 
= ,x,A,A*) is a set of five fields. The La- 

grangian density in Eq. (1) possesses a global U(l) sym- 
metry under the transformations 

^ e''^(t) , A e^'^A , (2) 

Adding sources terms J"(x) = {j,j*,s, S, S* ) to the La- 
grangian density, the generating functional W[J] of con- 
nected graphs is obtained from the path integral by [5] 



<s>,j]/h 



Z[J] = =M Jv<i> e"5' 

J] = + J d^x r{x) $„( 



(3) 



The generator of one-particle irreducible (1-PI) graphs, 
r[$], is obtained by a Legendre transformation from the 
classical currents J"(x) to the classical fields ^a{x) via 



r[$] = j d'^x r{x) $„(a;) - W[J] 



(4) 



The equations of motion and the inverse propagator are 



5r[$] 



Ja{x) , 



Q^^i^.A- (5) 



5^a{x) 5^^{x)5^'^{x') 

The path integration over the 4'{x) fields is done exactly 
and the integral over the fields x(x), A{x) and A*{x) is 
approximated by steepest descent. We obtain 

m ^lJJd^xd''x'4:ix)Gj[x,A]{x,x')M^') 



i^x{ 
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\A\' 
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— Tr 
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[HG-'[x,A]{x,x)]]} 



Here, G^j,^ represents the {1, 2} sector of G^p, he. 

Gj[x,A]{x,x') (6) 

ihdt-^v^ + x'r 
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where we introduced the notation x' = v^X ~ M- The 
finite-temperature partition function for the system is 
given by the same path integral (3) in the Matsub- 
ara imaginary time (t i— >■ —ihr) formahsm [5]. In the 
imaginary-time formahsm and for constant values of x' 
and A, the Green function Gablx^ depends only 

on X — x' . In Fourier space, we have 



V — UJ^t ] 



(7) 



where w„ ~ 2T:n/T are the Matsubara frequencies in 
fcs = 1 units. From Eq. (6), we find 



Cfe + X' + i^n 



A 



Gn(k,n) = G;2(k,n) = 
Gi2(k,n) = G2i(k,n) = ^ 2' 
with efc — and the dispersion relation 

The effective potential at finite temperature is[6] 

A(f>*^ (x' + A^)' 



(8) 



v,sm = x'\4>\' 



(27r)3 
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-Tlnd-e"'^"''"^ 



(9) 



(10) 



and the particle density is given by p = — dVcs/diJ = 
ix' + m)/(2A). The last term in Eq. (10) corresponds 
to a term Tr lnG~^ in the LOAF effective action [2]. 
Minimizing the effective potential, VcS, with respect to 
(j)* gives the minimum condition 



X>-A</)* =0. 



(11) 



Because of the gauge freedom, we can choose to be real 
at the minimum which means that for the broken sym- 
metry case ((^ 7^ 0), A = x' is real at the minimum, and 
leads to the Bogoliubov spectrum at weak coupling [2]. 
Minimizing the effective potential with respect to x' and 
A leads to the self-consistent equations: 
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= Po + A 



d^k l + 2n{u}k/T) 
(2^ 2^k 



P = Pa 



d^k ek + x' 



(2^)3 2ujk 



'l + 2n{ujk/T)] 



(12) 
(13) 



Here po — |0p and n{x) = {e^ — 1)~^ is the Bose distri- 
bution. Eq. (10) is regularized as shown in Ref. 4. 

The phase diagram in the LOAF approximation is de- 
picted in Fig. 1 as a function of the strength of the inter- 
particle interaction, which is characterized by the dimen- 
sionless parameter p^^^a, where ao is the s-wave scat- 
tering length. We notice the presence of three distinct 
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FIG. 1. (Color online) The LOAF phase diagram. 



regions, corresponding to the values of the three LOAF 
parameters, the usual (atom) BEC condensate density, 
Po, and the normal and anomalous auxiliary fields, x' 
and A: (I) < T < T^, with po 7^ and x' = ^ > 0, 

(II) Tc < T < T*, with Po = and x' > ^ > 0, and 

(III) T > T*, with Po = ^ = and x' > 0. Here, the 
critical temperature Tc corresponds to the emergence of 
the atom BEC condensate, whereas the temperature T* 
is related to the onset of superfluidity in the system and 
the emergence of a diatom condensate, A. The anoma- 
lous auxiliary field A represents a second order parameter 
in the LOAF theory. Region I in the phase diagram fea- 
tures both an atom BEC condensate and a superfluid 
state. In region II the atom BEC condensate is no longer 
present, but part of the system is still in a superfluid 
state, whereas in region III the entire system is in the 
normal state. In the noninteracting limit, and T* are 
the same. As the interaction strength increases, the tem- 
perature range for which the superfluid is present in the 
absence of the atom BEC condensates expands. LOAF 
predicts ^ 20% temperature range T* — Tc relative to 
Tc for a dimensionless parameter value, 

regions I and II fields carrying U{1) charge are nonzero, 
which leads to spontaneous breaking of the U{1) charge 
and the existence of Goldstone modes. These Goldstone 
modes are essential for the existence of superfluidity ac- 
cording to the Josephson relationship [7] . 

The superfluid state in the Bose gas is discussed us- 
ing Landau's phenomcnological two-fluid model [8]. We 
deflne the normal-state density 



Pn 



1 

3^ 



dk k^ ej. 



dn{ujk) 
duik 



(14) 



where uik is the quasi-particle energy given by the LOAF 
dispersion relation (9). Correspondingly, the superfluid 
density is defined as ps — p — Pn- In Fig. 2, we illustrate 
the temperature dependence of the atom BEC conden- 
sate, Po, and the superfluid density, p^, relative to the 
system density, p, for an interaction strength p^/^a = 0.4. 
As advertised, the onset of superfluidity in the system oc- 
curs at T*, and this temperature is different from the Tc, 
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the emergence temperature of the atom BEC condensate. 

If we promote the U(l) symmetry to a gauge symme- 
try by setting A gA{x)/h where g is the U(l) charge, 
in order to preserve the symmetry we introduce a gauge 
field Wij_{x) = ( Wo{a::), W(a:) ) analogous to the weak in- 
teraction vector gauge field. Then, the boson fields, </>, 
and the anomalous auxiliary field. A, carry one and two 
units of U(l) charge, respectively. The LOAF approxi- 
mation predicts that the superfluid state is accompanied 
by a Meissner effect in the presence of a weak vector 
potential, W(x). Following the standard derivation of 
the supercurrent in the Bardeen-Cooper-Schrieffer (BCS) 
theory (see e.g. Sees. 51-52 in Ref. 9) and adapting it to 
the case of a Bose gas of particles, one obtains the super- 
current in the low-momentum limit, as 



we obtain 



j.(q) = -p.— W(q) 
mc 



(15) 



Unlike in the BCS case, in LOAF the superfluid den- 
sity ps in Eq. (15) is the same as that obtained in the 
Landau two-fluid model described above. Characteristic 
signatures of superfluidity in Bose atom gases, such as 
dissipationless flow [10] and formation of quantized vor- 
tices in rotating gases [11], have already been observed. 
Superfluidity can also be measured directly by analyzing 
the "scissors" mode [12] with the atomic cloud oscillating 
with respect to the symmetry axis of the confining poten- 
tial. In nuclear physics this mode was observed in heavy 
deformed nuclei [13]. Recently, Cooper and Hadzibabic 
proposed using a vector potential generated by optical 
beams with nonzero angular momenta to simulate uni- 
form rotation of the atomic gas [14]. The induced change 
in angular momentum can be measured spectroscopically 
and provides a determination of the superfiuid fraction. 

The Goldstone theorem [15] states that when a con- 
tinuous symmetry, such as U(l), is spontaneously bro- 
ken, then, necessarily, new massless scalar states appear 
in the excitation spectrum related to the order param- 
eter. In LOAF, the composite-field Goldstone theorem 
for the auxiliary field A gives rise to a massless scalar 
when Tc < T < T*. To derive the Goldstone theorem, 
we consider Noether's theorem for the U(l) transforma- 
tion (2) on the fields <&. We calculate the change in the 
Lagrangian density, C = C + Ja^", under the infinitesi- 
mal change, $"(a;) ^'^{x)+S^°'{x), with (see e.g. Sec. 
19 in Ref. 16) 



S<^'^ix) = -eg%<^^{x), 



(16) 



where we 



introduced the U(l) charge matrix, 
g"/3 = g diag(l,-l,0,2,-2), for the fields 
= {0, 4i*,X: ^1 ^*} corresponding to a = 1 • • • 5. The 
complex fields (i>(x) and A{x) carry U(l) charges g and 
2(7, respectively. The real field x is U(l) neutral. From 



(5/:' 



sc 



dt 



+ V-jo(2;) - (iA) Ja5"/3*^ (17) 



where 



gh 



(18) 



2mi 



[(l)*{x)V(j3{x) - (\)(x)V<i)*{x)' 



Eq. (17) is a classical result, representing the U(l) current 
conservation in the absence of external sources. We mul- 
tiply (17) by exp[ — 5'[$, J]//i], divide by and integrate 
over the fields $ to derive the quantum version. Next, 
we substitute JaW> = Sr[^]/S^°'{x) from Eq. (5) in the 
resulting equation, where F is the Legendre transform of 
InZ and generates one-particle irreducible graphs. We 
integrate over d^a; and discard the surface terms to ob- 
tain the functional equation [16] 



S<^'r{x) ' 



(19) 



Differentiating (19) with respect to $"(3;) gives 



d X 



<52F[$] 



= / dV g-hx,x')g\<^P{x') = - J^(x) . 



For constant fields and in the absence of sources, Eq. (20) 
gives A4ap^^ — 0, where Aiap is a 4x4 matrix with 
indices {1, 2, 4, 5} and 

Map = I dV g-^{x,x')g^p = Qa^ima^'p • (21) 



The Goldstone theorem corresponds to det [ ] = and 
implies the presence of a pole in the propagator at zero- 
energy and zero-momentum transfer [16]. This result is 
preserved order-by-order in the auxiliary-field (AF) loop 
expansion, because the Ward identities are preserved 
order-by-order. In the LOAF approximation, the (j) prop- 
agator is leading order, but the mixing between the and 
A propagators only arises at next-to-leading order (NLO) 
and is proportional to 4>. Hence, det[ ] is given by the 
product of the determinants of the {1,2} and {4,5} di- 
agonal blocks, respectively, plus NLO corrections. The 
latter are dropped in LOAF. Then, the zero-determinant 
condition is satisfied if the determinant of either diago- 
nal block vanishes. Furthermore, at the minimum of the 
effective potential, the fields can be taken to be real, and 
the two determinant conditions are equivalent with 

fe' (0,0) -gr/(O,O)]0- 0, (22) 

[^4^(0,0) -^5:1(0,0)] A- 0. (23) 

Nontrivial solutions of Eqs. (22) and (23) require the 
brackets to vanish. These are the broken-symmetry Ward 
identities for the </> and A propagators, repectively. 
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FIG. 2. (Color online) (top) Comparison of the atom BEC 
condensate density, po, and the superfluid density, ps, for 



ao 



0.4. (bottom) Comparison of A and pa- Here 



Aq = A{T = 0). 



The Goldstone theorem for the atom BEC condensate 
corresponds to the case ^ 0. Hence, Eq. (22) is consis- 
tent with the minimum condition (11), which is satisfied 
in region I of the LOAF phase diagram: for T < Tc with 
X' 



A. In LOAF we find that in general 

^11(0,0) -gigl (0,0) 



(24) 



1 rl 



A 



d^fc l + 2n(cjfe/r) 



(27r)3 
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Here we used 

^m(o,o) = ^ 

^4:5(0,0) = ^ 



2 J (27r) 
d^fc 



(27r) 



/d fc v-^ ~ 
j^;^ 2Z <^ii(k, n) G22(k, n) 

3 ^ G'i2(k,7i)Gi2(k, n) , 



with the Green functions Gab given in Eq. (8). In region 
II {Tc<T < T*) where po is zero but A 7^ 0, Eq. (24) is 
identically zero and we find a composite-field Goldstone 
theorem, corresponding to a zero energy and momentum 
excitation of the gas. 

Next, we connect the superfluidity and the Higgs effect. 
From the supercurrent (15) and the Maxwell equation, 
V X (V X W) =}s, we find the London equation [17] 



W = 



(25) 



where \l ~ \J m(? j {^.-k p gg^^ is the London penetra- 
tion depth for the W field. To relate ps to the 
fundamental quantities of our theory, we note that 
the Tr{lnG~^} term in the action leads to nonlocal 
temperature-dependent n — (/) vertices that are poly- 
gons in the Green function, G. Following a similar ap- 
proach in BCS theory [18], we write a low-momentum 



Ginzburg-Landau effective theory for the 2(/-charge com- 
posite field, A. For simplicity, we use the relativis- 
tic generalization of our theory to make contact with 
the Anderson-Higgs mechanism [3]. By symmetry argu- 
ments, the relativistic effective field theory for the field A 
in the presence of a U(l) gauge field has the form 



C^{D^Ar{D'^A)~XA{\A\'~Ay 



with the covariant derivative, D,, — d, 



1 p pfii- 



(26) 



and 



Fn^ = daWu - d„Wn. Here, Aa 



^ - .^^ -r 2igWf„ 
= Xa{T) is the value of 
the four-point field interaction at zero momentum trans- 
fer and temperature T, and A = A(T) is obtained by 
solving Eqs. (12) and (13) for T < T*. Note that when 
g ^ 0, this is the usual effective Lagrangian for charged 
scalars which exhibits the Goldstone theorem. Using the 
gauge freedom, we expand the field about A as 

A^A + ^{Ai + iA2) , (Ai) = {A2) ^0, (27) 

which leads to the Lagrangian of the form 

C = '^idf,Aif-lXAA^Al (28) 
- I F^.F^" + ^{d^A2 + 2V2gAW^'Y +■■■ . 



Introducing a new field: W'^ 
the Lagrangian (28) becomes 



+ d^A2l{2^gA), 



c = \{d^A,f -\\aA'aI-\ f'^.f"^-- 

+ {2gAf W'^W'^" + ■■■ . 



(29) 



where F'^^ 



d^W'^. Hence, the field Ai has a 



composite-field Higgs mass, ~ XaA^, whereas the 
effective mass of the gauge field W'^ is — {2gA)^. 
The latter is identified as 



(29 A)' 



Ps 2 



(30) 



This implies that A^ is a measure of the superfluid den- 
sity, Ps- In Fig. 2 we show the temperature dependence 
of A'^ closely resembles that of the superfluid density. 

We note that one could also introduce U{1) gauge flelds 
into the BCS theory of dilute fermonic atom gases. In 
the fermion case, the gauge flelds couple to the diatom 
gap field. A, in much the same way as the gauge field 
couples to the composite- (diatom-) field A here. Then 
one would have derived the fact that the superfluid den- 
sity in BCS theory would be directly proportional to A^, 
as shown for Fermi systems in a more formal way from 
the Josephson relationship [7] by Taylor [19]. 

Equation (30) is consistent with the Josephson rela- 
tionship [7]. This fact is apparent in region II, where 
the derivation of the Josephson relationship parallels the 
corresponding calculation performed by Taylor [19] for 
Fermi systems. In contrast with LOAF, where the field 
A is composed of bosons, in BCS theory the gap, A, is 
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made up of fermions, causing a sign change in the loop 
contribution to the inverse propagator. Following Tay- 
lor's methodology [19], one obtains the superfluid density 
for a superfluid moving with velocity v, as 



Ps 



1 
V 



1 



v=0 



lim = , (31) 



where ^ly = —TlnZy is the grand canonical thermody- 
namic potential obtained by performing a twist opera- 
tion [20] on the A field, = e*^™^-^"'^, and is the 
A* A component of the composite field propagator. Just 
like in the BCS picture, in region II, ps ^ requires zero 
modes in the composite-field propagators. Apart from 
a wave function renormalization, Eqs. (30) and (31) are 
the same. We have shown recently [21] that the Joseph- 
son relation (31) holds also in region I, where (j) and 

To summarize, in this paper we studied the physical 
properties of a new phase predicted by the LOAF approx- 
imation to the dilute Bose gas theory. This phase resides 
in the phase diagram for temperatures between the criti- 
cal temperature Tc, where the BEG appears first, and the 
temperature T* corresponding to the onset of superfluid- 
ity. We proved that this phase supports composite-field 
(diatom) Goldstone states accompanied by zero-energy 
and zero-momentum excitations in the system. The pres- 
ence of superfluidity is tied to the existence of a U(l) 
charge-two diatom condensate in the system. We also 
showed that the mass of the gauge particle generated 
through the Anderson-Higgs mechanism is related to the 
superfluid density via the Meissner effect. We find that 
the square of the anomalous density is a good approxi- 
mation for the superfluid density. 
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